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We describe a method for the calculation of accurate energy eigenvalues and expectation values of 
observables of separable quantum-mechanical models. We discuss the application of the approach 
to one-dimensional anharmonic oscillators with symmetric potential-energy functions. 



I. INTRODUCTION 

Some time ago we developed a new method for the calculation of eigenvalues and 
eigenfunctions of the Schrodinger equation for separable problems. We first applied it to 
bound statesi'^j^i'^i^i^i^ but later we found that it also applied to resonances or metastable 
states^'^i'^'^'^, as well as to some strange complex eigenvalues^. 

Although the method only applies to separable quantum-mechanical problems it exhibits 
many interesting features that may be suitable for the discussion at a graduate, or advance 
undergraduate level. The teaching of the approach may be greatly facilitated by available 
computer algebra systems (CAS) that enable one to carry out the necessary calculation 
without much effort or special programming ability. 

In this paper we introduce the method and apply it to the one-dimensional Schrodinger 
equation with a symmetric potential. In order to make this article more comprehensive, in 
Sec. [TTlwe discuss the standard power-series treatment of the harmonic oscillator. In Sec. IIIII 
we develop the method from the known solutions of the harmonic oscillator. In Sec. IIVI we 
generalize it and apply it to some anharmonic oscillators with symmetric potentials. In 
Sec. |V] we discuss the phenomenon of multiplicity of roots that is one of the particular 
features of the approach. In Sec. I VI I we address the occurrence of meaningful spurious roots 
which is another peculiar feature of the method. In Sec. IVIII we consider the eigenfunctions 
and the calculation of expectation values. In Sec. we apply the method to a simple double- 
well oscillator. In Sec. [IX] we draw some conclusions on the teaching and utility of the 
procedure. 

II. THE HARMONIC OSCILLATOR 

In order to facilitate the introduction of the new approach and the understanding of most 
of its special features we first discuss an standard approach to the harmonic oscillator. For 
simplicity we assume that we have rewritten the Hamiltonian operator in dimensionless form 
as if = p^ + V{x) that in the coordinate representation reads H = —(P/dx^ + V{x). Thus, 
the Schrodinger equation for the harmonic oscillator reads 

^"(x) + {E- x'^)i){x) = 0, ^{x -^ ±00) = (1) 



The logarithmic derivative y{x) = ip'{x)/ilj{x) satisfies the Riccati equation y'{x) + y{xy — 
x'^ + E = 0. For large values of |x| we find that y(x) ~ ±x and ip^x) ~ e^^ /^. We appreciate 
that there are two solutions that we may call convergent ipd^) and divergent ipdi^) with 
asymptotic behaviours ^d^) ~ e~^ /^ and ipd{x) ~ e^ ^^ at large values of |x|. The most 
general solution is a linear combination of both: iIj{x) = Ci{E)iIjc{x) + C2{E)ipfi{x). The 
allowed values of the energy E are the roots Eq, Ei, ... of C2{E) = 0. Such eigenvalues En 
are consistent with a square integrable function /^ |^p(ix < cxo given in this case by ipdx)- 
A textbook method for the solution of the Schrodinger equation for the harmonic oscillator 
is based on a power-series expansion of the solution that we write asr^ 

^(x) = e-^'/2 J2 cjx^^+^ (2) 

3=0 

where s determines the behaviour oiiplx) at origin. The coefficients Cj satisfy the recurrence 
relation (2j + s + l){2j + s + 2)9+1 + [E - 1 - 2{2j + s)]cj = 0, and cq 7^ provided that 
s = or s = 1 that give rise to the even or odd states, respectively. 

One can prove that for arbitrary values of E the series in equation ([2]) behaves asymp- 
totically as e^ and the divergent contribution to the eigenfunction iIj{x) would dominate 
at large |x|. The eigenfunction will not be square integrable unless the otherwise infi- 
nite series in equation ([2]) becomes a polynomial, which occurs when E = En = 2n + 1, 
n = 21/ + s = 0,1, . . ., u = 0,1, . . .. When E = En then c^+i = and Cj = for all j > i^+l- 
The eigenf unctions then read 

V>„(X) = NnX'Pn{x)e-^"^^ Pn = JZc,,nX^'^' 

j=0 

'^■+''" ^ (2j + . + l)(2j + . + 2)'^ = °'^'---'^ ^^> 

The polynomials x^Pn{x) are proportional to the well known Hermite onesii, and the explicit 
form of the normalization factor Nn is unnecessary for our aims. 

The brief discussion given above has the sole purpose of calling the reader's attention 
on the close relationship between the correct asymptotic behaviour of the solution of the 
Schrodinger equation at |a;| ^ 00 and the allowed values of the energy, as well as on the fact 
that s/x — ^^(x) /'?/'„ (x) is a rational function, analytic at origin, and can be expanded in a 
Taylor series: 

^ YnK-^) _ X^n\X) — rn\-^J_ _ V^ r 2j+l /^\ 



X i)n{x) Pn{x) 



j=0 



III. THE RICCATI FADE METHOD FOR THE HARMONIC OSCILLATOR 

Suppose that we do not know the allowed values of the energy of the harmonic oscillator 
and consider the regularized logarithmic derivative of the eigenfunction 

that satisfies the Riccati equation 

f [x) + fw +x -E = (6) 

X x^ 

where the last term vanishes because s = or s = 1 as discussed earlier. In order to derive 

this expression we simply differentiate f{x) with respect to x and use equations ([T]) and ([5]) 

to get rid of ip" {x) / il){x) and ip' {x) / il){x) . According to what was said above we know that 

we can expand f{x) in a Taylor series about the origin 

00 
/(x)=x^/,■(E)^^ z = x^ (7) 

j=0 

Its substitution into the Riccati equation ([6]) leads to a recurrence relation for the coefficients 
fj (see below) that enables us to obtain as many of them as necessary. The first few 
coefficients of this expansion when s = are 

k = E 

E'-l 
h — 



f'. 
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2(^2 - 1)E 

15 
{E^-l)il7E^-5) 



f' = ' §15 («) 

Their general form appears to be fj = {E"^ — l)Qj{E) for j > 0, where Qj{E) is a polynomial 
function of the energy. 

Taking into account that the regularized logarithmic derivative of the eigenfunction is an 
exact rational function for the harmonic oscillator when E = En-, as shown by equation (jlj), 
we look for a solution of the form f{x) = x[M/N]{x^), where 

[M/iV](.) = lfi-1- = J: f,{E)z^+0{z'^^^^') (9) 



j:Ub,^' 



j=Q 



Since we can arbitrarily choose bo = 1 we are left with M + N + 1 coefficients of the rational 
function and the unknown energy as independently adjustable parameters; therefore we can 
try to obtain M + N + 2 exact coefficients of the Taylor series for f{x) as explicitly indicated 
in equation i^. When M > N we easily derive the following equations: 

v[iin(j,N) 

J2 hfj-k = ttj, j = 0, 1, . . . , M 

fc=0 

N 

Y.^kf,^k = 0, j = M+l,M + 2,...,M + iV + l 

fc=0 



(10) 



We can view the second set as a system of A^+1 homogeneous equations with A^+1 unknowns 
hp4, 6iv_i, . . . , &o- Therefore, there will be a nontrivial solution if i? is a root of 



H%{E) 



Jm-N+I fM-N+2 ■ ■ ■ /m+1 
fM-N+2 fA'I-N+3 ■ ■ ■ fM+2 







(11) 



Im+I fM+2 ■ ■ ■ Im+N+I 

where d = M — N = Q,l, . . . and D = A^ + 1 = 2, 3, . . . is the dimension of the Hankel 
determinant Hf){E). The ffist Hankel determinants for s = are: 



Hl{E) 
Hl{E) 
Hl{E) 
HUE) 



4725 
(^2 - 1)2(^2 -25)(E2 + 3) 

297675 

(^2 -1)3(^2- 25)2(^2 _ 81) 

46414974375 
4(^2 - 1)3(^2 - 25)2(^2 - 81)E 



(12) 



896041080309375 

The general form appears to be Hf) = {E^ - 1)^{E^ - 25)^"^ ...[E^ - {AD - 3f]G%{E) 
where Gj^{E) is a polynomial function of E. We appreciate that 

• As D increases more exact eigenvalues appear as roots of Hf)[E) = 

• The multiplicity of each root increases with D 

• In addition to the "physical" roots En = 2n + 1 there are spurious ones; for example: 
— (2n + 1) (later on we will discuss the occurrence of such negative roots) 

• We have obtained the allowed energies without taking into account the asymptotic 
behaviour of the eigenfunctions at large |x| explicitly 



IV. ANHARMONIC OSCILLATORS 

In this section we generalize the main results derived above for the harmonic oscillator 
and develop the method for the Schrodinger equation 

ij"{x) + [E- V{x)]-^{x) = (13) 

where V{x) is analytic at origin and symmetric about this point: V{—x) = V{x). Without 
loss of generality we assume that V{0) = 0. Since we can expand V{x) in a Taylor series 

oo 

V{x)=Y,VjX^^ (14) 

i=i 

then we can apply the method exactly as indicated above for the harmonic oscillator. The 
regularized logarithmic derivative of the eigenfunction f{x) = s/x — ip'{x)/ip{x) satisfies the 
Riccati equation 

/'(a;) + M^ _ f^^f ^ ^(^) _ ^ ^ (15) 

If we expand /(x) in odd-power series as in equation (jTj) we easily calculate the coefficients 
fj by means of the recurrence relation 

'" = 2n + L + 1 (g f'f"-'-' + ^''"'' - '''■) . " = 1- 2. ■ ■ ■ 

and then construct the rational approximation (j9]) that leads to the Hankel determinant 
(ITTj) . We do not expect to obtain exact eigenvalues for the general anharmonic oscillator 
as we have already done for the harmonic one, but at least we do expect that the roots 
of the Hankel determinants will enable us to estimate the eigenvalues of the Schrodinger 
equation flT3|) . We call this approach Riccati-Pade method (RPM) because it is based on 
the Riccati equation for the modified logarithmic derivative of the eigenfunction, and a 
rational approximation or Fade approximant i^^. 
The simplest symmetric anharmonic potential is 

V{x) =x^ + Xx'^ (17) 

This model has been widely studied because the perturbation expansion 

oo 

E{\) = Y,E^^^X^ (18) 

j=0 



is divergenip^ii^ii^ii^. A straightforward argument based on scaling the coordinate of the 
Schrodinger equation shows that the energy can also be expanded as 

oo 

Ei\) = \'/''f^e''^^\-''/^ (19) 

j=0 

with a nonzero convergence radiu9^>i^. 

From H2{E) we obtain (for concreteness we restrict to even states s = 0) 

E^ - 27E^ + 162E'^X + 51E'^ - 162E\ - 189A^ - 25 = (20) 

Fig. [1] shows that this simple expression yields reasonable results for the ground state of 
the anharmonic oscillator (TT7|) for all A values. The roots E{X) of (!20|) can be expanded 
in A-power series ( ITSl) and also in A^^'^-power series ( IT9l) . To prove the latter statement, 
simply substitute A^^^iy for E in equation (l20ll and notice that the roots of the resulting 
equation 

W^ - 27W*X-^^^ + 162W^ + blW'^X'*/^ - 1Q2WX'^/'^ - 25A"2 - 189 = (21) 

are functions of A^^/'^. Besides, when A — S' cx) we obtain 

W^ + 162W^ - 189 = (22) 

One of its roots W ~ 1.05 is quite close to the exact leading coefficient e*^°^ (calculated 
below) and explains why the simple expression ( l20l) yields satisfactory results for all values 
of A. 

From Hl{E) we obtain 

E^ - 2AE^ - 558E^X - 30E^ + UOiE^X + (eOGA^ + 128)^^ - 846^A - 882A2 - 75 = (23) 

It has been proved that roots oi H^{E) and H\){E) are lower and upper bounds, respectively, 
to the eigenvalues of the Schrodinger equation with the potential (TTTj) ^'^. Fig. [1] also shows 
the approximate ground-state eigenvalue given by Eq. ( l23l) and we appreciate that the 
upper and lower bounds are quite close for all values of A. 

In order to test the rate of convergence of the RPM we consider the most difficult case 
A -^ oo, sometimes called the strong-coupling limit. Since liniA^oo A^^/'^-E'(A) = e'-"-* is an 
eigenvalue oi H = p^ + x'^ we consider the potential 

V{x) = x^ (24) 



in what follows. Let E^^''^ be a sequence of roots of Hfj[E) that converges towards the 
energy eigenvalue En of the anharmonic oscillator with the potential ( l24l) . As said above, 
i5]j^'°l and E!^'^^ converge towards En from below and above, respectively, as D increases. 



p[D,l] p[D,0] 



that 



giving progressively tighter lower and upper boundsi>S. Fig. [2] shows log 
is a measure of the rate of convergence of the bounds towards the ground-state eigenvalue 
as D increases. From straightforward linear regression we estimate that Eq — Eq ^ 
77.31e~^'^^^, D = 2,3, . . .. Not many approaches exhibit such remarkable exponential rate 
of convergence. We find that E^^^'^^ = 1.0603620904841828996 is exact up to the last digit 
and that the accuracy of the RPM results increases with D according to the exponential 
rate of convergence just given. 

The number of zeros of ipni^) increases with n and for that reason we need determinants of 
greater dimension to take into account the increasing oscillation of the excited states (think 
of the rational approximation to f{x)). For the first excited even state n = 2 we obtain 
(8.02 X 10^) e-^-50^, D = 3,4, ... and Ef^'°^ = 7.4556979379867383922 is 



^[D,l] _ ^[D,0] 



exact up to the last digit. 



V. MULTIPLICITY OF ROOTS 



In the case of the harmonic oscillator we saw that the multiplicity of each root increases 
with the determinant dimension. The counterpart of this unusual feature of the RPM in 
the case of problems that are not exactly solvable is the occurrence of many roots in the 



p[D,0] p[oo,0] 



for the anharmonic 



neighbourhood of the chosen eigenvalue. Fig. [3] shows log 
oscillator with the potential f l2^ . where we have arbitrarily chosen Eq' = Eq for this nu- 
merical experiment. We appreciate that there are several parallel sequences of roots converg- 
ing towards the ground-state eigenvalue. A straight line log Eq — Eq' ^ 1.95(1 — D) 
marks the optimal sequence that was chosen above as a lower bound. The strategy for 
identifying the optimal sequence from the other ones consists simply in selecting the root of 
Hf){E) = closest to the one chosen previously from Hf)_^{E) = 0. 

For the first even excited state we find log E2 — E2' ~ 4.85 — 1.95D which is parallel 
to the ground-state one. 



VI. SPURIOUS ROOTS 

The RPM may also yield roots that are not related to the eigenvalues of the chosen model. 
To understand this interesting feature of the approach consider the change of variables x = 75 
that transforms H = —d?/dx^ + V{x) into H = 7"^ [—(P/dq^ + 7^V^(75')]. Therefore, some 
roots of the Hankel determinant for H will also give us an approximation to 7"^ times the 
eigenvalues of —d^/dg^ + 7^^(75) provided that 7^^(70') is a "reasonable" potential-energy 
functioni^. For example, in the case of the harmonic oscillator 7^V^(75') = 7^5'^ and the 
choice 7^ = —1 explains the occurrence of the negative roots — {2n + 1) discussed above. 

It is interesting to consider the modified Posch-Teller potential^ 



A(A-i : 

cosh(x) 



t'(-) = -^5ii (25) 



with exact eigenvalues En{MPT) = —{X — 1 — n)^. The change of coordinates discussed 
above with 7 = i yields —1 times the Hamiltonian operator with the Posch-Teller potential^ 

Vil) = ^^ (26) 

cos(gj^ 

and eigenvalues En{PT) = (A + 2n)^. We thus expect sequences of roots of the Hankel 
determinant converging towards En{MPT) and —En{PTY-. 

Table [T] shows them for n = and d = 0. This example is greatly revealing be- 
cause the asymptotic behaviour of the bound-state eigenf unctions of the first model fl25|) 
is tp{x —>■ ±00) = 0, whereas for the second one (|26ll the Dirichlet boundary conditions 
ilj{±7r/2) = are determined by the poles of the potential-energy function fl2^ at ±7r/2. 
We appreciate that the RPM yields solutions to two models with quite different boundary 
conditions simultaneously. This is partly a consequence of the fact that the RPM does not 
take the boundary conditions explicitly into account as discussed earlier in Sec. IIIII 

VII. EIGENFUNCTIONS 

Once we have a sufficiently accurate eigenvalue we obtain the coefficients fj and the 
rational approximation ([9]) to f{x). The approximate eigenfunction given by 

^(x) = Ara:^e-r/(3/)'^y (27) 



will be accurate in a neighbourhood of the origin but in general will not satisfy the boundary 
conditions at infinity. In principle, one can use the accurate energy and coefficients to derive 
an improved expression for the eigenfunction that behaves satisfactorily at infinity. However, 
if one needs to calculate the expectation value of an operator A it is preferable to apply the 
RPM to H + (3 A and calculate the slope dE/df3 at /3 = 0. According to the Hellmann- 
Feynman theorem we have^^ 

f(«-(i>t8) (28) 

Straigtforward differentiation of Hf){E,P) = with respect of /3 leads to dH'^/dp = 
{dHi/d(3)E + {dH'}^/dE)p{dE/dl3) = from which it follows that 

dE ^ idHi/d/3)E . . 

df3 (dHi/dE)^ ^ ^ 

Fig. m shows the convergence of the RPM values for (x^) for the ground state of the anhar- 
monic oscillator with potential 1^^ estimated by equations (128|) and (!29|) at /9 = 0. In this 
case we simply apply the RPM to the anharmonic oscillator with potential-energy func- 
tion V{x) = x^ + (3x^ and obtain (x^) = 0.3620226487886768452. It is worth noticing that 
(x^) = e*^^^ is the second-leading coefficient of the series flT^ . 

VIII. DOUBLE WELL OSCILLATOR 

The RPM developed above is also suitable for the calculation of the eigenvalues of the 
Schrodinger equation with symmetric double-well potentials. For concreteness and simplic- 
ity we consider the anharmonic oscillator with the potential-energy function 

V{x) = x* + I3x^ (30) 



that exhibits two wells when /3 < 0. They are located at x^ = ±J—P/2 and their depth is 
V{x^) = -/5V4. 

As in the case of the anharmonic oscillator discussed above we find that the RPM provides 
lower and upper bounds: -&„ '°^ < En < -£/„ '^ • We first consider the implicit analytic 
expressions given by H2{E,(3) = 0, and Hl{E,(3) = 0: 

^6 _ 27^4^ + iQ2E^ + 51E^p^ - 162Ep - 25/5^ - 189 = 
E^ - 2AE'^f3 - 558E^ - 30E^/3^ + UOAE^/3 
+ (666 + 128/3^)^2 - 846^/3^ - 882/3 - 75/3^ = (31) 

10 



respectively. Fig. [5] and [6] shows the roots of these equations as functions of p. We appreciate 
that one branch of each equation agrees with the accurate resuhs also plotted in those figures. 

Fig. Elshows log |eJ^'^^ - eJ^'°^ I for D < 20 and p = -1, -5, -10, -15, as well as straight 
lines provided by straightforward linear regression of the data for each value of (3. We 
appreciate that the rate of convergence decreases slightly and the RPM sequences start 
at greater D as f3 becomes more negative. However, in all those cases the convergence of 
the RPM is remarkable. We obtain similar results for the first excited state. The rate of 
convergence of the RPM bounds is considerably greater to that of the moment method^^. 
However, the latter method is more general and can be applied to a wider variety of problems. 

Table [IT] shows the first two eigenvalues obtained from Hankel determinants with d = 
0, d = 1, and D < 20. The upper and lower bounds agree to 20 digits for all those 
values of (3 except for (3 = -15. In this case we estimate -50.8413872843819547 < ^o < 
-50.8413872843819543 and -50.8413872841870053 < ^i < -50.8413872841870051. 

It is worth mentioning that present approach is based on a Taylor expansion of the 
solution about the top of the barrier, and it is well known that it is more convenient to 
expand about one of the minima. However, in the latter case we would have a nonsymmetric 
potential-energy function that we do not discuss here. 

IX. CONCLUSIONS 

In this paper we have presented an approach for the accurate calculation of eigenval- 
ues and eigenfunctions of separable quantum-mechanical problems. In this introductory 
communication we have restricted to one-dimensional models with symmetric (even-parity) 
potential-energy functions. The generalization to other cases is straightforward and will be 
discussed in subsequent articles. 

This presentation appears rather mathematical at first sight, but one must keep in mind 
that the Schrodinger equation for anharmonic oscillators have many applications in the 
study of, for example, oscillatory phenomena^i. Therefore, the RPM may be of practical 
utility in such studies, particularly because its remarkable rate of convergence enables one 
to obtain accurate results easily. We would not say that the RPM is the best choice for all 
possible applications, but it is worth taking into consideration that it does not require the 
calculation of matrix elements which may be difficult for some potentials and only requires 



11 



their Taylor expansions. Besides, in some cases, like the anharmonic oscillator (fT7|) . the RPM 
yields simple implicit equations that are reasonably accurate for all values of the potential 
parameters. 

But the main purpose of this paper is to discuss an approximate method that does not 
share many features with those presented in standard textbooks on quantum chemistry and 
quantum mechanics^i^i. It is those peculiar features of the RPM that in our opinion could 
make it most interesting for students with some knowledge in mathematics and physics. In 
particular, the RPM is suitable for training students in the use of CAS that facilitate most 
of the required algebra and analysis. 
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TABLE I: Ground-state eigenvalues of the models in equations (j25|) and (|26 



D 



^odlS 



^odM 



2 -3.9988835968549022343 -8.9983715148790601222 

3 -3.9999985313068943510 -8.9999994018534834235 

4 -3.9999999984549959290 -8.9999999999381908008 

5 -3.9999999999984803299 -8.9999999999999974948 

6 -3.9999999999999985472 -9.0000000000000000000 

7 -3.9999999999999999986 -9.0000000000000000000 

8 -4.0000000000000000000 -9.0000000000000000000 



TABLE II: Ground- and first-excited-state eigenvalues of the double-well oscillator (|30 
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0.65765300518071512306 


2.8345362021193042147 
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-3.4101427612398294753 


-3.2506753622892359802 
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-20.633576702947799150 


-20.633546884404911079 


-15 


-50.8413872843819543 


-50.8413872841870051 
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FIG. 1: Ground-state energy of the anharmonic oscillator (I17p . The solid and dashed curves are 
the analytical lower ()20p and upper ()23p bounds, respectively, and the circles mark more accurate 
results 
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FIG. 2: Convergence of the upper and lower bounds to the ground-state energy of the anharnionic 
oscillator ([24 




FIG. 3: Sequences of roots convergmg towards the ground-state energy of the anharnionic oscillator 
dMD as log 



p[D,0] p[«3,0] 
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FIG. 4: Convergence of the upper and lower bounds to (x^) for the ground state of the anharnionic 
oscillator with potential ([2] 
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FIG. 5: Roots E{(3) oi H2{E,j3) = (points) and accurate results (filled circles) for the ground state 
of the double-well oscillator ([50}) 
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FIG. 6: Roots E{P) oi H2{E , (3) = (points) and accurate results (filled circles) for the ground state 
of the double-well oscillator (1301) 




FIG. 7: Logarithmic convergence of the lower and upper bounds log \Eq — Ek \ for the ground 



state of the double-well oscillator (pOj) for /3 
(circles), /3 = — 15 (filled circles) 
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